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ABSTRACT 

We have surveyed a field covering 9.0 degrees 2 within the NOAO Deep Wide-Field Survey region in 
Bootes with the Multiband Imaging Photometer on the Spitzer Space Telescope (SST) to a limiting 
24 /mi flux density of 0.3 mJy. Thirty one sources from this survey with F24 Mm > 0.75 mJy which 
are optically very faint (R > 24.5 mag) have been observed with the low-resolution modules of the 
Infrared Spectrograph on SST. Redshifts derived primarily from strong silicate absorption features 
are reported here for 17 of these sources; 10 of these are optically invisible (R > 26 mag), with no 
counterpart in Bw, R, or /. The observed redshifts for 16 sources are 1.7 < z < 2.8. These represent 
a newly discovered population of highly obscured sources at high redshift with extreme infrared to 
optical ratios. Using IRS spectra of local galaxies as templates, we find that a majority of the sources 
have mid-infrared spectral shapes most similar to ultraluminous infrared galaxies powered primarily 
by AGN. Assuming the same templates also apply at longer wavelengths, bolometric luminosities 
exceed 10 13 L Q . 

Subject headings: dust, extinction — galaxies: high-redshift - infrared: galaxies — galaxies: starburst 
galaxies: AGN 



1. INTRODUCTION 

Opening the infrared wavelength regime for discov- 
ery is a primary objective of the final Great Observa - 
tory, the Spitzer Space Telescope l|Werner et alj [2004'). 
It was anticipated that new categories of sources would 
be revealed in various surveys. It was also anticipated 
that some of these sources would be very faint optically, 
as dust obscuration is probable in infrared-luminous 
sources, so the Infrared Spectrograph on Spitzer (IRS 7 - 
Houck, et al. 2004) was designed to include low resolu- 
tion modules which could determine redshifts for dusty 
sources. High-redshift, optically faint but infrared-loud 
;alaxies were first h inted at by the deepest IRAS surveys 
Houc ket alJll984tl . and partially uncovered in deep sur- 
veys with ISO (Dole et al. 2001, Elbaz et al. 2002). The 
sensitivity and wavelength coverage of the IRS enables us 
to measure the redshifts and power sources within such 
galaxies even if they are too faint for optical spectroscopy. 

In order to search for such optically-faint infrared 
sources, we have surveyed the Bootes field of the NOAO 
Deep Wide-Field Survey (NDWFS, Uannuzi fc Devi 
(1999)) with t he Multiband Im aging Photometer for 
Spitzer fMIPS: lRieke et alJ ((20041) ) to produce a catalog 
of mid- infrared sources in 9.0 deg 2 to flux density lim- 
its at 24 /jm of 0.3 mJy. This field was chosen because 
the deep and well calibrated optical imagery in B\y, R, 
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and I bands makes possible the identification of infrared 
sources with very faint optical counterparts and confident 
selection of infrared sources lacking optical counterparts 
to very deep limits. Our overall objective for the cur- 
rent study was to select the infrared sources which are 
faintest optically but bright enough for feasible redshift 
determination with the IRS, taking a limit of 0.75 mJy 
for these initial observations. 

Details of the survey will be reported elsewhere. We 
used the MIPS and optical surveys to select and inspect 
all sources with 24 /im flux densities above 0.75 mJy and 
with I > 24 mag. Of the 4273 MIPS 24 /im sources 
brighter than 0.75 mJy, 114 met this optical magnitude 
criterion. Of these, 27 sources have no optical iden- 
tification in any band, meaning they are fainter than 
about 26.5 mag in Bw, 26 mag in i?, and 25.5 mag in 
I (all Vega magnitudes). Based on comparison of posi- 
tions between the many identified optical counterparts 
and infrared sources, the formal 5 er positional criterion 
for non-identification of 24 /jm sources with flux densities 
above 0.75 mJy is a distance > 2" from the nearest op- 
tical source. We also rejected infrared sources close to a 
bright galaxy or within clusters of faint galaxies with the 
presumption that redshifts of the infrared sources can 
be determined by association with the optical sources. 
As a result of these additional criteria, our initial target 
list contained 17 sources with no optical identifications, 
and all of these have been observed with the IRS. We 
also observed an additional 13 sources having an optical 
counterpart but fainter than 24.5 mag in R. Finally, one 
source was included (#13 in Table 1) because it is rel- 
atively bright at 24 /im but has an extreme IR/optical 
ratio (2.3 mJy at 24 /im and i?=23.8mag). In this pa- 
per, we report redshifts, physical properties and lumi- 
nosities for 17 of the 31 sources which show sufficiently 
strong spectral features that redshift measurements can 
be made with the IRS. 
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2. OBSERVATIONS AND DATA REDUCTION 

Observations and results for the sources discussed in 
this paper are summarized in Table 1, giving the name 
(with full coordinates), the 24 /im flux density, optical 
magnitudes of the MIPS 24 /im source if detected, in- 
tegration times for the IRS observations, redshift deter- 
minations, and other source characteristics. The MIPS 
data were obtained with an effective integration time of 
~90sec per sk y pixel and re duced with the MIPS Data 
Analysis Tool {Gordon et afll2004|) . Point source extrac- 
tion was performed using an empirical point spread func- 
tion (PSF) constructed from the brightest objects found 
in the 24 /im image and subsequently fitted to all the 
sources detected in the data. We allowed for a multiple 
match fitting to deal with blended cases. We finally de- 
rived the flux of each object using the scaled fitted PSF 
after applying a slight correction to account for the finite 
size of the modeled point spread function. 

Spectroscopic observations were made with the IRS 
Short Low module in order 1 only (SL1) and with the 
Long Low module in orders 1 and 2 (LL1 and LL2), 
described in iHouck et al.l {2004). These give low reso- 
lution spectral coverage from ^8 /im to ~35 /im. Sources 
were placed on the slits by offsetting from nearby 2MASS 
stars. Observed images were processed with version 10.5 
of the SSC pipeline. Background subtraction and ex- 
traction of source spectra alo ng the slit was done with 
the SMART analysis package {Higdon et al.ll2004bl) . All 
spectra discussed in this paper are shown in Figures 1 
and 2. 

3. DETERMINING REDSHIFTS OF SOURCES 

We determine redshifts by using bright, local "tem- 
plate" sources for which we have full spectral coverage 
with the IRS. In almost all redshift fit depends on 

the presence of a strong decrease in the continuous spec- 
trum beginnin g around 9 /im in the source rest frame (for 
example, see lHigdon et al.l {2004alU . This drop is identi- 
fied with the strong silicate absorption feature centered 
at 9.7 /jm ; this feature is common in ultraluminous in- 
frared galaxies (ULIRGs) at low redshift (e.g. Armus 
et al. 2004). An extreme case of this type of absorp- 
tion is seen i n the IRS sp ectrum of the ULIRG IRAS 
F00183-7111 {Spoon et alJl2004aj) . Strong spectral emis- 
sion features, such as the 7.7/xm and 8.6 /*m PAH fea- 
tures, can be present without the presence of deep silicate 
absorption - e.g., in the unobsc ured starburst nucleus 
of NGC7714 (Brandl et al. 2004). The well known low- 
redshift ULIRGs Arp220and Mrk231 are intermediate 
cases. Arp 220 has deep silicate ab sorption and strong 
PAH emission {Spoon et al.ll2004bl) . while Mrk231 has 
weak silicate absorption, but little or no PAH emission. 
To cover the range from pure absorption to pure emis- 
sion, we have fit all spectra with the IRS templates of 
NGC7714, F00183-7111, Arp 220, and Mrk231. (The 
IRS spectra of Arp 220 and Mrk 231 will be discussed in 
more detail in future papers.) A formal \ 2 routine gives 
the optimal redshift fit for each template. Once such fits 
were formally determined, the four were compared to es- 
timate which one gave overall the best fit to the observed 
spectrum. That is the result listed in Table 1 and shown 
in Figure 1. 

Figure 2 illustrates the effects of alternative templates 
for our brightest 24 /im source, #9 in Table 1, showing 



that the redshift uncertainty for a given source primar- 
ily depends on source S/N and strength of the features 
rather than choice of template. For the sources in Table 
1, this uncertainty typically is ± 0.2 in z; the formal un- 
certainty derived from the \ 2 fit f° r the template adopted 
is given for each source in Table 1. 

The 17 sources out of our 31 observed spectra with 
measured redshifts derived from template fits are all 
shown in Figure 1. The 14 remaining sources do not 
have sufficiently strong features to list a spectroscopic 
redshift. An example of one of these sources is also shown 
in Figure 1 to illustrate the obvious difference between 
spectra with and without strong mid-infrared features. 
Of the 17 sources, 15 have unambiguous redshifts based 
on the x 2 fits shown, and all of these have 1.7 < z < 2.8. 
The two sources for which we have less confidence in 
the redshift include the only source, #11, with assigned 
z < 1. The redshift of 0.70 arises because all three of 
the templates containing absorption give a better fit at 
this redshift than at the higher redshifts characterizing 
all of the other sources; however, the spectrum beyond 
30 /xm for this source does not have sufficient S/N to 
show the drop in flux density which would be present for 
a strong absorption feature at high redshift. Source #2 
is the only source for which the template fits produce an 
ambiguous redshift based on % 2 . Comparably good fits 
arise from Arp 220 at z = 0.57 or NGC 7714 at z = 1.86. 
The NGC 7714 fit is adopted because of our presumption 
that a dusty source is more likely to be optically faint at 
high redshift, but if this fit is correct, this is the only 
source with a pure PAH emission spectrum. 

4. DISCUSSION 

These results clearly demonstrate the presence of a 
population of optically-obscured infrared sources at red- 
shifts 1.7 < z < 2.8. There are significant selection effects 
in our result that constrain the redshift range in which we 
can identify this obscured population. We have deliber- 
ately selected against optically visible sources, which gen- 
erally forces sources to significant redshift, probably z > 
1, because the rest-frame UV is most affected by internal 
obscuration. For 1.0 < z <1.6, any deep 9.7 /un silicate 
absorption feature in the spectrum would significantly 
affect survey detections at 24 /im, so absorbed objects 
would selectively fall out of the 24 /im sample. For z > 
3.1, the continuum drop indicating the absorption feature 
would move out of the IRS spectral coverage. It is not 
surprising, therefore, that the redshifts we have identified 
fall in the range 1.6 < z < 3.1. Within this range, our 
initial results represent only a lower limit for the extent 
of this obscured population. This is primarily because 
of our restriction to sources brighter than 0.75 mJy at 
24 /im, but also because many sufficiently bright infrared 
sources with faint optical counterparts (the 84 remain- 
ing sources having an optical identification but with I 
> 24 mag) were excluded from IRS observation in the 
hope that redshifts for these very faint sources might be 
determinable with ground-based spectroscopy. 

Even though the sources in our sample are faint, 
their high redshifts require luminosities that place them 
among the class of extreme or "hyper-luminous" infrared 
galaxies. However, these sources have been chosen be- 
cause they are extremely faint in the optical, imply- 
ing much larger intrinsic extinction than in a typical 
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ULIRG. For example, an object having the infrared SED 
of Arp 220 and sufficiently luminous to have a 24 /im flux 
density of 1 mJy when at z = 2 would have / <23 mag, 
making it too bright in the optical to be included in our 
sample. (A source of the same luminosity as Arp 220 if 
at z — 2 would be much too faint to be included in our 
infrared sample, having a 24 /im observed flux density of 
only 0.03 mJy.) Thus, while our sources have IRS spec- 
tra that can be well fit with known ULIRG templates in 
the rest-frame mid-infrared, they are all much fainter op- 
tically, and therefore apparently more heavily extincted, 
than low-redshift ULIRGs. What is the nature of this ob- 
scured population? Is the mid-infrared luminosity of the 
hot dust produced primarily by an AGN or by a starburst 
in our sample galaxies? If IR AS F00183-7111 is po wered 
by an AGN, as suggested bv lSpoon et al.1 (|200 4b"l . then 
13 out of 17 (76%) of the sample are best fit with an AGN 
mid-IR template (F00183-7111 or Mrk 231). Four of the 
galaxies are best fit with a starburst template (Arp 220 or 
NGC7714). Although the fraction of light coming from 
an AGN or from a starburst in any ULIRG is uncertain, 
these template fits suggest that most of the sources we 
have found are AGN dominated - i.e. steeply rising spec- 
tra with weak PAH emission features and varying levels 
of silicate absorption. 

Radio data can also aid in understanding the nature of 
these sources. A well defined correlation exists between 
infrared and radio flux densities for syste ms with no AGN 
contri bution to the radio flux density. lAnpleton et al.l 
(2004) show that this relation as described by the pa- 
rameter q = log[f„(24 /im)/f l/ (20 cm)] has q = 0.8 ± 
0.3. The Bootes survey area has been mapped as p art of 
a Westerbork 20 cm survey ()de Vries et al. 1 120021. Al- 
though the radio beam (13" x 27") is much larger than 
even the MIPS 24 /im beam, we attempt source iden- 
tifications by considering any 24 /jm source that has a 
quoted Westerbork source within 5" to be an identifica- 
tion with this radio source. With this criterion, five of the 
infrared sources have radio counterparts. For those that 
do not, we take the upper limit to the 20 cm flux density 
to be a 5<7 value of 0.20 mJy. The q values are listed in 
Table 1. Except for source #13, all sources have values 
or limits for q very close to those expected for starbursts 
or radio-quiet AGN. Source #13 does have the excess ra- 
dio emission which would characterize a radio-loud AGN, 
and this source is best fit with the IRAS F00183-7111 
template, confirming at least for this case the validity of 
the AGN interpretation. 

To understand the contribution of this obscured pop- 
ulation to the luminosity density of the universe at high 
redshift, estimates of bolometric luminosities are needed. 
Our only way to estimate these is to assume that mid- 
infrared luminosities scale to bolometric luminosities the 
same as in the template sources which fit the mid- 
infrared spectra. Using this scaling, all sources in Ta- 
ble 1 for which we have z > 1.7 have implied infrared 
luminosities of 6 x 10 12 - 6 x 10 13 L Q , with an average 
value of 2.3 x 10 13 L©. The largest existing samples of 
dusty sources with comparable bolometric luminosities 
and redshifts are t he SCUBA or MAMBO submillimeter- 
selected galaxies l)Chanman et al.l 12003). Detection at 



24 /tm of such sources brighter than a few mJy at submil- 
limeter wavelengths is common, but nearly all of these 
galaxies are fainter at 24 /im by a typica l factor of 3- 
5 than the sources r eporte d here (e.g. Ilyison et al. I 
lEffll. lEgami et all ll200l . iFraver et all l|2004fl . and 
Charmandaris et al. ( 2004)). This implies a higher ratio 
of mid-infrared to far-infrared luminosity in our sources 
than in the sub-millimeter sources. This indicates that 
our new sample seems to represent sources that have 
hotter dust than typical submillimeter-selected galaxies; 
this could be an additional component of hot dust or 
an overall averaged dust temperature which is greater. 
Either would be consistent with having AGN as the pri- 
mary luminosity source instead of the st arbursts thought 
to pow er the sub-millimeter galaxies. iChaoman et alJ 
(2004) utilized optically determined redshifts to iden- 
tify a population of optically faint radio galaxies (R > 
23.5 mag) at z ~ 2 which are not detected as submillime- 
ter sources, and half of this sample has optical indicators 
of an AGN. Alt hough none of the sour ces were observed 
in the infrared, IChaoman et all l)2004|) suggest from the 
absence of submillimeter detections that this sample is 
also characterized by hot dust, so it may represent a sim- 
ilar population to our sample. 

5. SUMMARY 

The NOAO Deep Wide-Field Survey area in Bootes 
has been surveyed at 24 /im with the MIPS instrument 
on SST, and 31 optically faint or invisible sources have 
been observed with the IRS instrument. Our first re- 
sults indicate that sources chosen solely on the basis 
of extreme 24 /im to / band flux density ratios, ^f„(24 
/zm) / 'ui v {I) > 100, are heavily selected in favor of red- 
shifts exceeding 1.5 and that the majority of such sources 
have redshifts that are measurable with IRS based on 
a strong silicate absorption feature in the mid-infrared. 
Sixteen of the 31 objects selected for the initial IRS ob- 
servations are found to have spectroscopic redshifts 1.7 
< z < 2.8, based on fits of redshifted IRS spectral tem- 
plates for local ULIRGS, and 9 of these have no optical 
counterparts. The template spectra and radio survey 
data suggest that this population is dominated by AGN- 
powered, radio-quiet ULIRGs. If bolometric luminosities 
scale with mid-infrared luminosities as in the templates, 
these newly discovered sources have bolometric luminosi- 
ties between 6 x 1O 12 L and 6 x 10 13 L Q . 



We thank many others on the IRS, MIPS, and Bootes 
Survey teams for contributing to this effort, especially 
including D. Barry, D. Devost, P. Hall, G. Sloan, H. 
Dole, P. Appleton, and C. Bian. This work is based in 
part on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided 
by NASA through Contract Number 1257184 issued by 
JPL/Caltech. This work used data provided by the ND- 
WFS, which is supported by NOAO. NOAO is operated 
by AURA, Inc., under an agreement with the NSF. 



REFERENCES 



Appleton, P. N., ct al. 2004, ApJS, 154, 147 



Armus, L., et al. 2004, ApJS, 154, 178 



4 



Houck et al. 



Brandl, B. R., et al. 2004, ApJS, 154, 188 

Chapman, S.C., Blain, A.W., Ivison, R.J. and Smail, I.R. 2003, 
Nature, 422, 695 

Chapman, S.C., Smail, I., Blain, A.W. and Ivison, R.J., 2004, ApJ, 
616, 71 

Charmandaris, V. et al. 2004, ApJS, 154, 142 i 
de Vries, W.H. et al. 2002, AJ, 123, 1784 
Dole, H., et al. 2001, A&A, 372, 364 
Egami, E. et al. 2004, ApJS, 154, 130 

Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, H., Franceschini, 

A., Fadda, D., & Chary, R. R. 2002, A&A, 384, 848 
Frayer, D.T. et al. 2004, ApJS, 154, 137 
Gordon, K., et al., 2004, PASP, submitted 
Higdon, S. J. U., et al 2004, ApJS, 154, 174 
Higdon, S. J. U., et al 2004, PASP, 116, 975 



Houck, J. R., et al. 1984, ApJ, 278, L63 
Houck, J. R., et al. 2004, ApJS, 154, 18 
Ivision, R.J. et al. 2004, ApJS, 154, 124 

Jannuzi, B. T. & Dey, A., 1999, in "Photometric Redshifts and the 
Detection of High Rcdshift Galaxies", ASP Conference Series, 
Vol. 191, Edited by R. Wcymann, L. Storrie-Lombardi, M. 
Sawicki, and R. Brunner. ISBN: 158381-017-X, p. Ill 

Rieke, G. H., et al. 2004 ApJS154, 24 

Spoon, H. W. W., et al. 2004, ApJS, 154, 184 

Spoon, H. W. W., Moorwood, A. F. M., Lutz, D., Tielens, 
A. G. G. M., Siebenmorgen, R., & Kcanc, J. V. 2004, A&A, 
414, 873 

Werner, M. W., et al. 2004, ApJS, 154, 1 



High Redshift Obscured Sources 



5 



TABLE 1 

Redshifts and Source Characteristics 



Source 




Name a 




f v 24/im 
(mjy) 


Bw 
(mag) 


R 

(mag) 


I 

(mag) 


timc b 
(sec) 




z 




Template 


IR/opt c 


q 




T d 

(xl0 13 L 


1. 


SST24 


J142958.33+322615 


4 


1.24 


25.6 


25.1 


24.4 


1200 


2 


64+ 


25 


F00183-7111 


110 


>0 


8 


1.8 


2. 


SST24 


J142653.23+330220 


7 


0.89 






24.7 


1680 


1 


86+ 


07 


NGC 7714 


100 


>0 


8 


0.9 


3. 


SST24 


J143447.70+330230 


6 


1.81 








960 


1 


78+ 


30 


Mrk 231 


>270 





8 


0.9 


4. 


SST24 


J143523.99+330706 


8 


1.08 


26.7 


24.6 


23.5 


1440 


2 


59+ 


34 


F00183-7111 


40 





8 


1.5 


5. 


SST24 


J142804.12+332135 


2 


0.87 








1680 


2 


34+ 


28 


F00183-7111 


>130 


>0 


6 


0.9 


6. 


SST24 


J143358.00+332607 


1 


1.03 








1440 


1 


96+ 


34 


Mrk 231 


>150 


>0 


7 


0.6 


7. 


SST24 


J143251. 82+333536 


3 


0.78 


26.2 




24.3 


1920 


1 


78+ 


14 


Arp 220 


60 


>0 


8 


2.8 


8. 


SST24 


J143539.34+334159 


1 


2.65 


26.2 


25.4 


24.5 


720 


2 


62+ 


26 


F00183-7111 


250 


>1 


1 


3.8 


9. 


SST24 


J143001. 91+334538 


4 


3.83 








480 


2 


46+ 


20 


F00183-7111 


>570 


1 





4.6 


10. 


SST24 


J143520.75+340418 


2 


1.45 


25.5 


25.3 


24.3 


1200 


2 


08+ 


21 


Mrk 231 


110 


>0 


8 


1.0 


11. 


SST24 


J143242.51+342232 


1 


0.87 








1680 





70+ 


21 


Mrk 231 


>130 


>0 


6 


3.0 


12. 


SST24 


J142626.49+344731 


2 


1.12 








1440 


2 


13+ 


09 


Arp 220 


>170 





8 


5.7 


13. 


SST24 


J143644.22+350627 


4 


2.30 


24.7 


23.7 


23.4 


720 


1 


95+ 


17 


F00183-7111 


80 


-0 


3 


1.8 


14. 


SST24 


J142538.22+351855 


2 


0.79 








1920 


2 


26+ 


11 


Arp 220 


>120 


>0 


6 


5.3 


15. 


SST24 


J142645.71+351901 


1 


1.05 








1440 


1 


75+ 


21 


Mrk 231 


>160 


>0 


8 


0.5 


16. 


SST24 


J142924.83+353320 


3 


1.04 








1440 


2 


73+ 


19 


F00183-7111 


>150 


>0 


7 


1.6 


17. 


SST24 


J143504.12+354743 


2 


1.23 








1200 


2 


13+ 


17 


Mrk 231 


>180 


>0 


8 


0.9 



a SST24 source name derives from discovery with the MIPS 24^tm images; coordinates listed are J2000 24^tm positions with typical 3<t 
uncertainty of ± 1.2"; sources with an optical counterpart will also appear in NDWFS catalogs with prefix NDWFS and the optical source 
position; optical magnitudes are Mag-Auto from NDWFS Data Release Three, available at http://www.noao.edu/noao/noaodcep/ |. 

b Integration time for each order of Long Low spectrum; integration time in Short Low order 1 is 480 s in all cases. 

c IR/opt = iT„(24 /im)/i/f„(7). Limits are based on assumed / > 25. 

d The instrinsic infrared luminosity of each source, Li r .=L(8-1000/xm), is estimated by scaling the known L; r of the best template fit to 
the measured MIPS 24/mi flux density after redshifting the template SED to its measured z. 
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10 15 20 25 30 35 10 15 20 25 30 35 10 15 20 25 30 35 

Wavelength (u.m) Wavelength (u.m) Wavelength (^m) 



Fig. 1. — Observed spectra smoothed to approximate resolution of individual IRS orders (solid line), best template fit (dashed line), 
observed flux densities and rcdshift for all 17 sources in Table 1. Spectra truncated below 15 /Ltm had no detectable signal in SL1; spectrum 
#14 also had no signal in LL2. Final panel is a source not listed in Table 1 to illustrate a spectrum without sufficient features for a rcdshift 
measurement. 
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_g ,'\Z = 237- 


g Z = 2.46: 
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'tem^J^FO^^ 
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10 15 20 25 30 35 10 15 20 25 30 35 10 15 20 25 30 35 
Wavelength (um) Wavelength (|xm) Wavelength (urn) 



Fig. 2.— Spectra of source #9 showing alternative template fits with Arp220, absorbed ULIRG IRAS F000183-7111, and Mrk231. 
Range of z values illustrates uncertainty in z that arises from different template assumptions. 



